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Abstract 

A new solvable many-body problem of goldfish type is identified and used to revisit the connection 
among two different approaches to solvable dynamical systems. An isochronous variant of this model 
is identified and investigated. Alternative versions of these models are presented. The behavior of 
the alternative isochronous model near its equilibrium configurations is investigated, and a remarkable 
Diophantine result, as well as related Diophantine conjectures, are thereby obtained. 
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I. Introduction 

Recently a method has been introduced and exploited to identify new exactly solvable (namely solvable 
by purely algebraic operations, such as diagonalizing a matrix) many-body problems characterized by 
equations of motion of Newtonian type ("the acceleration of each particle is determined by the positions 
and velocities of all particles"), including in particular models of goldfish type (see, for instance, 0], and 
below). The main idea of this approach - hereafter referred to as the direct method - is to start from 
an explicitly solvable matrix evolution equation (possibly even quite a trivial one), and to then focus 
on the time evolution of the eigenvalues of this matrix. For an overview of this method (including an 
explanation of the terminology used herein), of the main results yielded by it so far, and the quotation 
of relevant references, we refer to the very recent paper 2 . 

Another method has been introduced some years ago to treat certain well-know solvable dynamical 
systems and to illuminate their connection with developments in theoretical particle physics. The main 
idea of this approach - hereafter referred to as the gauge theory method - is to start from a gauge 
invariant matrix evolution equation and to exploit the possibility that in one gauge this evolution be 
trivially simple hence solvable while in another gauge it be related to interesting evolutions, in particular 
to the equations of motion of Newtonian type of certain many-body problems. For an overview of this 
approach and the quotation of relevant references, we refer to the relatively recent paper £Q. An analysis 
of the gauge theory approach entailing a clarification of the relation of this method to the direct approach 
is already provided in the more recent paper published by one of us (EL) [TT] , 

In the present paper, in the context of revisiting this connection, we identify a new solvable many-body 
problem of goldfish type. This finding hinges on a result obtained many years ago by V. I. Inozemtsev [ID) . 
We also present the isochronous variant of this many-body problem, as well as alternative formulations 
of these two models, and by investigating the behavior of the alternative isochronous model in the 
neighborhood of its equilibrium configurations we identify certain remarkable Diophantine relations. 

The main new results obtained in this paper are reported in the following Section II. The hasty 
browser eager to see immediately the equations of motion of the new solvable many-body problems of 
goldfish type should jump to and for the isochronous variant to ©, and for the alternative versions 
of these models to ()14jl and 1)21(1 : a Diophantine finding and related conjectures are reported at the end 
of Section II. In Section III the solvable character of the new many-body problems of goldfish type is 
demonstrated, firstly via the direct method and then via the gauge theory method; the connection among 
these two approaches is thereby illuminated. In Section IV solvable dynamical systems are derived, which 
constitute nontrivial alternative reformulations of the many-body problems of goldfish type treated in 
Section III. In Section V the behavior of the alternative isochronous model in the neighborhood of its 
equilibrium configurations is investigated and remarkable Diophantine relations are thereby obtained. 
In Section VI possible future developments are mentioned. The Appendix contains some findings the 
insertion of which where they are first mentioned (see Remark 2.9 in Section II) would have been too 
distracting. 



II. Main results 

In this section we report the main new findings obtained in subsequent sections. 

The solvable A-body problem of goldfish type identified in this paper is characterized by the following 
equations of motion of Newtonian type: 



Notation: z n = z n (t) are the dependent variables, t is the independent variable ("time"), superim- 
posed dots denote time-differentiations, a 2 is an arbitrary constant (we use a 2 rather than a merely for 
notational convenience, see below), A is a positive integer (generally we assume A > 1), and indices such 
as n, m generally take all the values 1, 2, A unless otherwise mentioned. 

Remark 2.1. Trivially related models involving additional arbitrary constants could of course be 
obtained by rescaling the (dependent and independent) variables and by shifting by a constant amount 




(1) 
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the dependent variables; note incidentally that the first factor 2 in the right-hand side of Q could be 
changed by rescaling (we put it there for notational convenience, see below), while the second factor 2 
(that multiplying the sum) cannot of course be changed. □ 

Remark 2.2. Although for real a 2 and for real initial data z n (0), z n (0) the time evolution (for real 
time) of this many-body model entails that the dependent variables z n (t) are as well real, we generally 
assume the time evolution to take place in the complex z-plane (and generally allow the constant a 
to be as well complex); indeed such an evolution is much more interesting due to the possibility of the 
"particles" characterized by the complex coordinates z n (t) to go round each other and the related fact that 
initial data z n (0), i„(0) leading to particle collisions are then exceptional (they generally have vanishing 
dimensionality relative to generic initial data). (If attention is instead restricted to real motions, then 
the trivial change of dependent variables z n — > i y n with y n real might be expedient in order to deal with 
confined motions.) It is possible to reformulate these complex equations of motions as real (and even 
covariant, even rotation-invariant) equations of motion describing the motion of real point particles in 
the real (say, horizontal) plane, but we will not take space here to reformulate them in this manner, since 
the technique to do so is well-known (see for instance Ref. [5J)- 

The solvable character of these equations of motion is evidenced by the well-known fact ^H] that 
the N x N matrix evolution equation 

U = 2U(U 2 -a 2 ) (2) 

is itself solvable (in terms of appropriate sigma functions [TJ]]), together with the following 

Proposition 2.3. The solution of the initial-value problem for the equations of motion Qjj) is provided 
by the following prescription: the coordinates z n (t) are the N eigenvalues of the N x N matrix U (t) 
solution of 0) and determined by the following initial data: 

U nm (Q) = S nm z n (0) , (3a) 

Unm(0) = Snm [z„ 2 (0)-a 2 ] 

+ [z„(0) + z 2 (0) -a 2 } 1/2 [i m (0) + zl(0) -a 2 } 1 ' 2 . □ (3b) 

Note that the matrix U(0) is diagonal, while the matrix ?7(0) is the sum of a diagonal matrix and a 
dyadic matrix. 

Notation: here and hereafter 5 nm = 8 n>m is the Kronecker delta symbol, 5 nm = 1 if n = m, S nrn = 
if n m. 

To obtain the isochronous variant of this many-body problem one starts from the equations of motion 

C = 2<n + 2 J2 + C ^^' m + C "') , (4) 

which correspond to Q with a — and with the merely notational replacement of the dependent variables 
z n (t) with the dependent variables £„ (r) (and of course now the appended primes denote differentiations 
with respect to r). One can then apply the procedure usually referred to as "the trick" (see for instance 
0), i- e. (in this case) the following change of dependent and independent variables 

z n (t) = cxp(it)C n (r) , (5a) 



r = i [l-exp(it)] . (5b) 

This yields the equations of motion 



N [z n -iz n + z 2 \[ z m - i z m + z^ 



S n = 3iz n + 2Z n (l + z 2 ) + 2 V ^ L. (6) 



The solution of the initial-value problem is then obviously given by the solution (via Proposition 
2.3) of the problem (J2J and by the "trick" relations (J5J, that clearly also imply 

C„(0) = z~„(0) , C'(0)=5„(0)-i5„(0) . (7) 
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Equivalently, the solution of this model I© is clearly given by the following 

Proposition 2.4. The dependent variables z n (t) that solve the initial-value problem for the Newtonian 
N-body problem Jj5J) are the A eigenvalues of the A x A matrix U(t) evolving according to the solvable 
matrix evolution equation 

U -3iU -2U = 2C/ 3 (8) 
and being moreover characterized by the following initial data: 

Unm(0) = S nm Z n (0) , (9a) 



EWO) = S nm [z 2 (0)] 



z n (0)-iz n (0) + ~z 2 {0) 



1/2 



z m (0)-iz m {0) + z„(0) 



1/2 



□ (9b) 



Note that the matrix U(0) is diagonal, while the matrix U(0) is the sum of a diagonal matrix and a 
dyadic matrix. 

The solvable character of the matrix evolution equation (JSJ is implied by the " trick" formula 

U(t) = cxp(i t) U{t) , r = i[l- exp(i t)] (10) 

relating the A x A matrix U(t) evolving according to JSJ) to the A x A matrix U(t) evolving according 
to J2J with a — 0. 

Remark 2. 5. The solvable character |10| of the matrix evolution equation (J2J entails that all its solu- 
tions U(t) are meromorphic functions of the independent variable t. Hence (see (llOfl ) all the nonsingular 
solutions U (t) of the matrix evolution equation arc periodic with period 2 tt, 

U(t + 2n) = U{t) . (11) 

The singular solutions of (JSJ are exceptional, corresponding to a set of initial data having vanishing 

measure with respect to the set of generic initial data. □ 

As an immediate consequence of Proposition 2.4 and of this Remark 2.5 there holds the following 
Proposition 2.6. All the solutions of the many-body problem (except those exceptional ones that 

run into a collision of two or more particles, which correspond to nongeneric initial data) are completely 

periodic with a period which is a positive integer multiple p of 2 n: 

z n (t + 2pn) = z n (t) , p = l or 2 or.. .or N . (12) 

The positive integer p accounts for the possibility that the eigenvalues get exchanged among each other 
through the motion: it depends on the initial data, but it does not change for sufficiently small, if finite, 
changes of these data and it clearly is not larger than N . □ 

This proposition displays the isochronous character of the A^-body problem indeed it justifies 
considering it as one more instance of nonlinear harmonic oscillators [7]- 

There exists a, by now rather standard, technique to reformulate these type of A-body problems, by 
identifying the N "particle coordinates" z n (t) as the A zeros of a (monic) polynomial in z of degree A, 
and by then focusing on the corresponding time evolution of the iV coefficients c m (t) of this polynomial 
(see for instance jE] [H]) : 

N N 

i,(z,t)=Hlz-Z n (t)] = J2 C rn( t )z N ~ m > C = 1 ■ (13) 
n— 1 Tci— 

In Section IV we show how such a procedure is applicable in our case, and we thereby obtain the following 
alternative formulation of the A-body problem Q: 

Cm + 2 (m - 1) c m+ i - 2 ci c m + 2 (A + 1 - m) a 2 c m -\ 
+ (to + 2) (m - 3) c m+2 - 2 (to - 1) c\ c m+ i 
+2 [to (A + 2 - to) a 2 + ci - c 2 + 3 c 2 ] c m 

-2 (A + l - to) a 2 c 1 c m _i + (N + 2 - to) (A + 1 - to) a 4 c m _ 2 = , 
to = 1,...,A, c = 1 , c_i = c N+1 = c N+2 = . (14) 
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Remark 2. 7. The ODE of this system with m = is identically satisfied; the ODE with m = N + 1 
is also satisfied provided one sets ctv+3 = 0, and even the ODE with m = N + 2 is identically satisfied if 
one moreover sets cn+a — 0. □ 

Remark 2.8. A superficial look at this system of ODEs might suggest that it is a linear system of 
evolution equations for the quantities c m (t); but this is of course not the case, due to the presence of 
the quantities c\(t) and C2(t). Indeed the highly nonlinear character of this system is already evident 
by looking at the N = 2 case, in which case it yields the following (solvable!) fourth order ODE for 
f(t) = Cl (t): 

f"" f 2 ~ 2 /"' /' / 2 - 2 /"' / 3 - 2 (/") 2 / + 2 /" (/') 2 + 4 /" /' / 2 - 2 /" / 4 

-4 (/') 2 / 3 +4/'/ 5 +4 a 2 (/"/ 2 -2/'/ 3 ) =0 (15) 

(here for typographical convenience differentiations are denoted by appended primes rather than super- 
imposed dots). □ 

Remark 2.9. Two equilibrium (namely, time- independent) solutions of this system (|14H are provided 
by the formula 

c m = {±a) m (^j. (16) 

They are not, however, the only equilibrium configurations of this system. A technique to obtain all 
these configurations (including this one!) is described in the Appendix. □ 

As indicated above, see l)13[). the quantities c m (t) that evolve according to the system of ODEs (|14|l 
are just the coefficients of the monic polynomial ip(z,t) of degree N in z, the N zeros z n (t) of which 
evolve according to the equations of motion (JTJ. Hence (see Proposition 2.3) the solution of the system 
of ODEs l|14f> is given by the following 

Proposition 2.10. The dependent variables c m (t) that solve the initial-value problem for the system 
of nonlinear ODEs {1$ are the N coefficients of the polynomial t/j(z,t), see flty). which is itself given by 
the formula 

ip(z,t) = det[z-U{t)} , (17) 

where the N x N matrix U(t) evolves according to the solvable matrix evolution equation 0) and is 
moreover characterized by the initial data 0), with the initial values Z n (0), z n (0) related to the initial 
values c m (0),c m (0) by the formulas implied by I J.'j|). 

N N 

l[[z - z n (0)} = ]T c m (0)z N - m , c = 1 , (18a) 

71 — 1 771 — 

N N N 

-$>„(0) [] [z-z n (0)] = J2^)z N - m .B (18b) 

n=l m =l,m=£n m=l 

To obtain an alternative version of the isochronous A-body problem © we use the following version 
of the "trick": 

c m {t) = (~i) m exp(mit) 7 m (r) , (19a) 



T = i[l- exp(it)] . (19b) 

Here the quantities 7 m (r) are the dependent variables of the previous model, ljl4[) . with a = 0, up to 
the (purely notational) change consisting in calling the independent variable r (instead of t) and the 
dependent variables 7 m (instead of c m ), so that these variables satisfy the following system of ODEs: 

7™ + 2 (m - 1) 7™+i - 2 7i 7m + (m + 2) (m - 3) 7m+2 
-2 (m - 1) 7i7 m+ i + 2 ( 7 i - 7i 2 + 3 72 ) 7 m = 
m = 1,...,A, 7o = 1 , 7_! = 7w+1 = 7^+2 = , (20) 

where of course appended primes denote differentiations with respect to the independent variable r (which 
we allow to be complex, see H19bjl ), 
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Then clearly by applying the "trick" 1)19)1 to the system (|20Jl the following new system of nonlinear 
ODEs is obtained: 



c m + 2 (to — 1) i c m+ i — (2 ?7i + 1 + 2 c\) i c m 

- (to + 2) (to - 3) c m+2 + 2 (to - 1) (to + 1 + e"i) c m+i 

+ -to (m + 1) + 2 i ci - 2 (m - 1) ci + 2 &f - 6 c 2 

m = l, ...,iV, c = 1 , c_i = cjv+i = cjv+2 = . 



= 



(21) 



Remark 2.11. The prefactor {—i) m in (|19a|l is of course unessential, it has been introduced merely to 
give a marginally nicer look to this system (|21|l and to some other formulas, see below. With this version, 
(|19|l , of the " trick" the relation among the particle coordinates satisfying the equations of motion of the 
isochronous TV-body problem © and the quantities c m (t) satisfying this system of ODEs 1)21(1 reads now 



JV 



N—m 



CO 



(22) 



see (0, 1)190 and (|13[) . Note that we introduced here the (new) monic polynomial ip(z,t) having as its N 
zeros the N dependent variables z n (t) satisfying JBJ and as its N coefficients the N dependent variables 
c m (t) satisfying □ 

This model, 1)21)1 . is obviously just as solvable as the previous one, 1)14)1 . indeed the solution of its 
initial-value problem can be obtained from the solution of the corresponding problem for 1)14)1 via the 
formulas 1)19)1 that clearly imply the following relations among the initial data of the two models: 



5m(0) = H) m 7m (0) , 



(23a) 



c m (0) - to i Cm (0) = H) m 7m (0) . (23b) 

Equivalently, the solution of this model 1)21)1 is clearly given by the following 

Proposition 2.12. The dependent variables c m {t) that solve the initial-value problem for the system 
of nonlinear ODEs \21\) are the N coefficients of the polynomial ip(z,t), see which is itself given by 

the formula 



ip{z,t) = det z - U(t) 



(24) 



where the N x N matrix U(t) evolves according to the solvable matrix evolution equation (JHJ and is 
moreover characterized by the initial data with the initial values z n (0), z n (0) related to the initial 
values c m (0),c m (0) by the following formulas implied by 



N 



N 



n^-^(o)]=Ew 



,(o). 



n=l 



N 

E 



z n (0) -iz n (0) 



m=0 
N 



N 



n [z-5„(o)]=e« 

m=l,m^n m— 1 



Cq = 1 



.(0). 



□ 



(25a) 



(25b) 



As an immediate consequence of this Proposition 2.12 and of Remark 2.5 there holds the following 
Proposition 2.13. All the nonsingular solutions of the system of ODEs 121}) are completely periodic 
with period 2 n, 

c rn (t + 2Tr) =c m (t) , (26) 

while the singular solutions are exceptional, corresponding to a set of initial data having vanishing mea- 
sure with respect to the set of generic initial data. □ 

This proposition displays the isochronous character of the iV-body problem 1)21)1 . indeed it justifies 
considering it as one more instance of nonlinear harmonic oscillators [Jj. 
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Finally, in Section V we obtain all the equilibrium configurations of the isochronous systems (JHJ) and 
(|21|l and we study the behavior of the system of nonlinear harmonic oscillators (|21|l in the neighborhood 
of its equilibrium configurations. The interested reader will find these results in that section, but we 
advertise here the Diophantine findings arrived at via this study. 

Proposition 2.14. Let the two N x N matrices A and B be defined componentwise as follows: 

A nm = 2 (n - 1) <5„+i, m - (2 n + 1 + 2 ci) <5 n , m + 2c n 5i, m , (27a) 



B nm = (n + 2) (n - 3) 6 n +2, m — 2 (n — 1) (n + 1 + ci) S n +i.n 
+ [n (n + 1) + 2 (n - 1) ci - 2 + 6 c 2 ] <5„,™ 
+2 [- (n - 1) c n+ i + (n - 1 - 2ci) c„] <5i, m + 6c n S 2 ^ 



with the numbers c m defined as follows: 

for v = 



for v 
r in 



(-)' 



5om + $1 



if /i = 1 , 







c 


--\) 


\ m J 









if /i > 1 , 



for ^ = 3 



(-)" 



m 



— 3 

771 — 1 

4 



for i/ = 4 , c m = (-) m 



for u = b 



(-)' 



/i — 5 

771 — 5 



fe=0 
5 



E 

fc=0 



+ 14 

771 — fc 

/i — 5 

777 — fc 



/j — 3 

777 — 2 

5 
k 



14 



/i — 3 

777 — 3 



c arbitrary 



(27b) 
(28a) 

(28b) 
(28c) 
(28d) 
(28e) 



^4s indicated above the parameter v (the role of which here is mainly to distinguish 5 different cases) can 
take any one of the 5 values 0, 1, 3, 4, 5, Tw/ii/e the parameter fi can take any positive integer value in the 
range v < jj, < N. Let the 2 N numbers p^ be the eigenvalues of the generalized eigenvalue problem 



(p 2 + Ap + B) r = , 
(where r= (ri, rjv) denotes the corresponding eigenvector) implying 



N 

det (p 2 + Ap + B) = JJ 



7J - p. 



(+) 



(29a) 



(29b) 



Then the 2N numbers p„ are all integers. □ 

Notation: here and throughout the symbol (^) is the standard binomial coefficient, 

r(x + i) 



yj r(y + l)r(x-y + l) ' 



(30) 



We have verified with the help of symbolic programing languages (we used Maple and Mathematica) 
and for an ample sample of values of JV and of the other parameters the validity of this proposition (proven 
in Section V), and from these computer-aided checks we are led to formulate the following Diophantine 
conjectures. 

Conjecture 2.15. For v = 0, 1, 3, 4, 5 and p, integer in the range v < p, < N the eigenvalues of the 
generalized eigenvalue problem ]29a)) (with {27}) and f28p ) are given by the following formulas: 

fori/ = 0, det (p 2 + Ap + B) 
= 1 JJ [(p-») 0-n-l)]l |jj[(p + n) (p + n-5)}\ , (31a) 
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for v = 1 , det (p 2 + Ap + B) = (p + 1) (p - 4) 
J] [(p-n) (p-n-5)] m[(p + n) (p + n-7)]l 



(31b) 



for 1^ = 3, det (p 2 + Ap + 5) = (p+1) (p - 4) ■ 
f] [(P-n) (p-n + 5)] ^ [] ttP + n ) (p-n + p-t)] , 



, 71=1 



(31c) 



for i/ = 4 , det (p + Ap + -B) = (p+ 1) 



][[ [(p-n) (p-n+1)] 



>-4 



J] (p + n) 



JJ(p_„_l) 

_n=l 

[J (p + n + 1) 



(31d) 



for i/ = 5 
i f A" 



det (p 2 + Ap + 5) 



H [(p-n) (p-n-1)] ^ [] [(p + n) (p-n + p-4)] 



(31e) 



ffere we use i/ie standard convention according to which a product equals unity if the lower limit of the 
running index exceeds the upper limit. □ 

Remark 2.16. For p = v = the validity of this conjecture is certainly true, indeed trivially so (see 
below the Remark 5.4). H 

The Conjecture 2.15 only refers to integer values of the parameter p in the range v < p < N . But 
our computer-aided exploration also indicates the validity, for arbitrary values of the parameter p, of 
the following conjecture (which is only formulated below for sufficiently large values of N, to avoid less 
interesting complications). 

Conjecture 2.17. The generalized eigenvalue problem $29a\) (with and y28a\) ) features, for 

arbitrary p, the N — 1 eigenvalues 



2, 3, 4, 5 - p, 6 - p, . 



N - 



if v = or v = 5 and N > 5 , 



(32a) 



2, 3, 4, 4 - p, 5 - p, AT - 1 - p , if v = 4 and N > 5 
and the N — 4 eigenvalues 

-1, 4, 6, 8 - p,9 - p, N - p , ifv= \ andN>&, 



(32b) 
(32c) 



-1, 4, 6, 3-p,4-p,...,7V 



if v = Z and N > 8 . □ 



(32d) 



Remark 2.18. The Conjecture 2.17 - in contrast to the Conjecture 2.15 — does not provide the 
entire spectrum of the eigenvalue problem l|29al) . which of course features 2N eigenvalues. □ 

Some aspects of these conjectures are easy to prove. For instance Conjecture 2.15 can be proven 
by induction for all N > p if one assumes its validity for N = p. But complete proofs of them do not 
seem quite trivial (see Section VI). 



III. Two proofs of Proposition 2.3 

In this section (part of) the results reported in the preceding Section II are proven, firstly by the direct 
method, then by the gauge theory method. 
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A. Direct method 



The starting point of the direct method is the solvable N x N matrix evolution equation (J2J. We then 
introduce the eigenvalues of the matrix U (t) and the corresponding diagonalizing matrix R(t) via the 
formulas 

U(t) = R(t)Z(t) [RW 1 , (33a) 



Z (t) = diag [z n (t)] , (33b) 

with moreover 

R(0) = 1 . (33c) 



Here and below 1 is the N x N identity matrix. Note that the first two of these equations, Ij33a(l and 
(|33b(l . identify (consistently with Proposition 2.3) the N coordinates z n (t) as the N eigenvalues of the 
N x N matrix U(t), while the third, (|33c|l . is consistent via the first two with the assignment (I3all . 

It is then easy to see (for the derivation of these formulas see, if need be, for instance [2]) that, after 
introducing the N x N matrix M (t) via the assignment 



one gets 



entailing (see (|33cl 



M(t) = [R(t)} 1 R(t) , (34) 
U(t) = R(t) [z(t) + [M(t), Z(t)]} [R (t)}- 1 , (35a) 



(7(0) = Z(0) + [M (0), Z(0)] , (35b) 

as well as the following system of evolution ODEs for the coordinates z n (t) and for the matrix elements 
M nm (t) of the matrix M(t): 



N 



= 2 z n (z, 2 - a 2 ) + 2 J2 ( z « - z ™) M «™ M ™« > ( 36a ) 



-2 Zn _ Zm - M nn + M„ 



N 



r, 



+ 2J z _ z — ^f^, n^m. (36b) 

Note that the time evolutions of the diagonal elements M nn (t) of the matrix M(t) remain unrestricted: 
it is indeed clear from 133fl (implying that R(t) is defined only up to multiplication from the right by an 
arbitrary diagonal matrix D(t)) and from Ij34fl that these N functions of time can be chosen arbitrarily 
without affecting the eigenvalues of U(t), namely the coordinates z n (t). Indeed it is clear that the N x N 
matrix evolution equation ©, characterizing the time evolution of the N 2 matrix elements M nm (t), has 
now been turned into the system H3(jfl , characterizing the time evolution of the N coordinates z n (t) and 
the N (N — 1) off-diagonal elements M nm (t) (with n 7^ to) of the N x N matrix M{t). 

Clearly this system l|36|) is no less solvable than the original matrix evolution (0) , because its solution 
can be retrieved from the solution of |J3J) by purely algebraic operations (essentially, by diagonalizing an 
N x N matrix). 

But we are interested in obtaining an TV-body problem involving only the N "particle coordinates" 
z n (t), hence our next task is to eliminate the N (N — 1) "auxiliary quantities" M nm (t) (with n =/= m). 
To do this one must find (assuming it exists) an appropriate ansatz expressing the N (N — 1) auxiliary 
quantities M nm (t) (with n =/= to) in terms of the N particle coordinates z n (t), taking advantage if need 
be of the freedom to assign the N quantities M nn (t) at our convenience. 

An ansatz that works (in the sense of turning the N (N — 1) evolution equations (|36bJ) into identities) 

is 

JV 

M nn (t) = \ ^ , (37a) 
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M nm {t) = - — — — — ^ , n^m , (37b) 

[z n (t) - z m [t)\ 

with g an arbitrary constant. This leads however to an iV-body model the solvable character of which 
is already well known [5] J3| HH Q2]> hence we do not pursue this development here (we elaborate this 
point a little further in the following version of the proof) . 

Remark 3. 1 . Clearly insertion of this ansatz (|37(1 in the more general matrix evolution equation 

U=$(U) , (38a) 

with $(z) an arbitrary scalar function would also work (since this does not depend on the equations of 
motion satisfied by the coordinates z n , see below), and it would lead to the iV-body problem characterized 
by the Newtonian equations of motion 

N 2 

Z n = 4> (z n ) - 2 f Z ^3 ■ (38b) 

This was already noted, many years ago, by Veselov But it appears that, so far, the most general 
(up to trivial transformations) solvable N x N matrix evolution of type l|38ajl is just (J2J). □ 

Another ansatz that also does (as it were miraculously) work (namely, transform the evolution equa- 
tions l|36b(l into identities) reads as follows: 

M mm (t) = - {[ *" + *"- 02] [^ + ^-° 2 ]} 1/2 > (39 ) 

Zn z m 

Note that this ansatz, in contrast to the previous one, contains no arbitrary ("coupling") constant g. In 
this case the appropriate assignment for the diagonal elements M nn (t) is quite trivial: M nn (t) = 0, or 
equivalently (see l|36b|0 M nn (t) = fi(t), fi(t) being an arbitrary function of time (but independent of the 
index n). The truth of this assertion can be verified by a trivial if tedious calculation: note that the 
evolution equations (|36a|l must also be used in the process. 

And it is now clear that the insertion of this ansatz in 13(ial) yields (|TJ, while its insertion in (|35b|l 
yields the assignment (|3b|l . The proof of Proposition 2.3 is thus completed. 



B. Gauge theory approach 

Let us now prove again Proposition 2.3, but via the gauge theory method. Although this entails some 
repetitions we believe it is useful to go through this exercise in some detail, especially because we will 
now use a somewhat different language - and one purpose of this paper is precisely to clarify the relations 
among these two different approaches. Moreover this presentation provides some indication of the extent 
to which this kind of fishing expeditions are likely to yield new goldfishes, namely new interesting solvable 
models. 

Let us start by reviewing (but in a notation more conducive to a direct comparison with the preceding 
treatment) the gauge theory approach to a more standard model referred to in the literature as rCM pQ 
[TT] , characterized by the Newtonian equations of motion 

N 2 

z n =-z n -2 J2 J-^ 73 ■ ( 4 °) 

We then present a variant of this approach which leads us to the new solvable iV-body model of goldfish 
type presented above. 

We start with the following system of dynamical equations for the N x N matrices U = U(t), M = 
M(t), and W = W(t): 

U + [M, U]=W , (41a) 
W + [M, W] = $([/) , (41b) 
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with the square brackets indicating matrix commutators. As discussed below, to get the rCM model one 
should assign the function $(£/) as follows: 



$([/) = -U . (42) 

It is however convenient to leave this function $ unspecified for the moment; but we require that it 
contain no other matrix besides U, so that <£(£/) R = ^(Rr 1 U R) for any (invertible) matrix R. 

The equations 1)41(1 have a natural interpretation as gauge theory in + 1 dimensions: they are 
indeed of the form [D t ,U] — P, [D t ,P] — with D t — dt + M being the so-called covariant time 

derivative with M regarded as gauge field. In particular, they are invariant under the following gauge 
transformations, 

U -> U = R- 1 UR, W-+W = R~ 1 WR , (43a) 
M -> M = R- 1 MR + R- 1 R, (43b) 

where the matrix R = R(t) characterizing the gauge transformation is an arbitrarily time-dependent 
invertible matrix (the transformation rule for M follows from D t R — dt + R^ 1 R + R^ 1 M R). One 
can exploit this invariance to impose additional conditions. In particular for any solution U(t), M(t), and 
W(t) of H41|) . one can find a gauge transformation R(t) such that the gauge-transformed matrix M(t), 
see l(43b(> . vanishes, M(t) = 0. Indeed, this is implied by the fact that the linear first-order matrix ODE 



R + MR = 0, R(0) = 1 (44a) 
always has a (unique) solution, which can be written as 



R(t) = T exp 



- / dsM(s) 
Jo 



(44b) 



where the symbol T denotes time-ordering. This shows that we can impose the condition 

M(t) = (45) 

for all times t, without loss of generality. Note that this also implies that, if we impose l|45H . we will 
not loose any solution: indeed any solution to our gauge theory equations can be obtained by solving 
these equations with the condition (|45|l imposed, and performing a gauge transformation afterwards. In 
particle physics this latter condition is often called Weyl gauge. We will also use another gauge condition, 
namely 

U nm (t) = S nm z n (t) (46) 

for all times t, which we call diagonal Coulomb gauge. Note that we can impose this latter condition if the 
matrix U(t) is such that there exists an invertible matrix R(t) such that R~ 1 (t)U(t) R(t) is a diagonal 
matrix, and this is obviously true in the generic case when the matrix U(t) is non-degenerate. The cases 
when the matrix U(t) is degenerate correspond to particles in our A^-body system colliding, and then our 
solution breaks down (as it should): but this can only happen for exceptional - i. e., nongeneric - initial 
data, if we allow the particle coordinates z n (t) to move in the complex plane, as we generally do (see the 
Remark 2.2). 

The idea now is that, by imposing the gauge condition (|4T)l . we get the matrix equation U = $({/) 
which might be (chosen to be) exactly solvable, whereas by imposing the condition (|46J1 . we get a (hope- 
fully) interesting dynamical system for the variables z n (t). The latter dynamical system can then be 
solved explicitly as follows: one first determines the solution U(t) of the matrix equation obtained from 
(|41|l in the Weyl gauge and with the initial conditions 

(0), U nn (0) = i„(0). (47) 

Then the eigenvalues of U (t) give the solution of the dynamical system. Note that we can only assign 
the diagonal elements of f (0) since, as we will see, the off-diagonal elements of £7(0) are determined by 
another conditions which we have to add. 
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Indeed, to get an interesting dynamical system, we need to add one more gauge invariant equation to 
(|4*TJ) . In particular, to get the rCM model one has to add the equation 

[W, U] = J (48) 

which is often called Gauss law or momentum map. This latter equation is gauge invariant if the matrix 
J = J{t) introduced here transforms under gauge transformations as J — > J = R^ 1 JR. It turns indeed 
out pQ that if one makes the assignment 

Jnm(O) = 9 (1 - 8 nm ) (49) 

and chooses $(£/) as in (|42(l . then the coordinates z n (t) obey the equations of motion of the rCM model 
IP. 

To obtain goldfish type dynamical systems one must instead replace the Gauss' law condition l|48|l by 



Bnra B n ' m ' — B nm ' J3 n ' m (50a) 

with 

B = W + f(U), (50b) 

where f(x) is another function to be assigned later. To see that this condition is gauge invariant we 
note that it can be written asB®B — VB^B where <g> is the tensor product (so that (B <g> B) (u ® 
v) = (B u) <S> (Bv) where u and v are iV-vectors) and V is the permutation matrix defined as follows, 
V u®v = v<£)u: the gauge invariance of (|50a|l follows from the obvious fact that V commutes with R®R. 
A new finding (proven below) is then given by the following 
Proposition 3.2. The gauge theory equations and $50\) are consistent if 

f(x)=a + l3x + jx 2 , (51a) 



$(x) = {a + (3x + 1 x 2 ){[3 + 2 1 x) =f(x)f'(x) , (51b) 
for arbitrary constants a, (3, 7. Imposing the diagonal Coulomb gauge J^6| ) these equations imply 

z„ = $(z„)+ J2 % + f(z n )][z m + f(z m )] ^ (52) 

■ / Zn Z-ui 

and thus the solution of the initial-value problem for this dynamical system can be obtained by solving the 
corresponding gauge theory equations in the Weyl gauge, see \4-Q^j . More specifically: the solution of the 
initial-value problem for this dynamical system, \5ty. is given by the eigenvalues of the matrix equation 

U = <P(U) (53) 

with the initial conditions 

(0) (54) 

and 

U nn (0) = z n , (55a) 



U mn (0) - {[i„(0) + / [z n (0)}} 1/2 {[i m (0) + / [z m (0)]} 1/2 , n ± m . □ (55b) 

Remark 3.3. In the special case 7 = a model is obtained whose solvability was already known 
[S]. Our new model reported in Section II (and already derived by our other method in the preceding 
Subsection A) is obtained for (3 — 0, 7 = 1, and a = —a 2 . But the greater generality of the result as 
formulated in Proposition 3.2 is only apparent: if 7 ^ 0, one can always reduce this more general case 
to the special case with (3 = by the (rather trivial) transformations U — * U = U — 7^-l,t —> t = jt. And 
note that, for f(z) — z 2 — a 2 , the right-hand sides of (|55|l and i|3b|l coincide: the apparent differences 
are merely notational. □ 
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Let us end this section by outlining the proof of this result, whose analogy with that proven in the 
first part of this section is we trust evident enough not to require further elaboration. We firstly write 
out H41f) and Q5U|I by imposing the condition (|45|l . Then (|41a|l becomes 

0~nm Zn ~f~ M nm (z m Z n ^j ^W n m (^6) 

which for the diagonal elements (i.e. n = m) implies 

W nn = i„ . (57) 
For the off- diagonal elements we obtain the assignment 

M nm = - W nm , nytm, (58) 

Zn Z ra 



whereas the diagonal elements M nn remain unassigned. Then (|41b(l reads 



N 



W nm + W t™ ~ M t™) = 6 «™ $ ( Z «) ' ( 59 ) 



The diagonal elements of this equation give (via (|B^|l and JHZ|)) 

W n iWl n 

2^ Z = ®\ Z n), (60) 

while the off-diagonal elements imply the following important consistency conditions, 

W nm + (M n - M m ) W nm + Zn ~ Zm W nm 

N ( i l \ ( 61 ) 

- V W n i Wi m = , n ^ m. 



Z n - Zl 



In particular (|5Ua|l implies B nm B mn — B nn B mm - Inserting in this equation the assignment B nm = 
Wnm + S nm f ' ( z n) (see H50b|l ) and in particular B nn = z n + f(z n ) (sec l|57)l). we get 

W nm W mn + 8 + f(Zn)} f(z n ) = [i„ + f{Zn)] [Z m + f{z m )] ■ (62) 

Thus the solution of (|5U|I for n' = m and m' = n is 

1/2 

W nm = exp(</?„ - (p m ) {-6 nm [2i„ + f(z n )} f{z n ) + [Z„ + /( Zn)\ \Zm + f{z m )]Y" (63) 

where the functions f n (t) are arbitrary, and it is easy to check that this is a solution of 15U|) also for all 
the other values of n',m'. Inserting this in (|6L)[) we obtain 

z„ = $(z„) + 2 J2 [i " + /(2 " )] _ [i£ + /(Z£)] ■ (64) 

n i Zn Zl 

It remains to check the consistency relations 1)61(1. We note that, for £ ^ n,m and n ^ m, 16311 entails 
Wni Wim = W nm \h + f(ze)], hence (ED is implied by 

N 



£ [it + /(%)] ( — — ) , n ? m . (65) 



e=l,£^n,m 



Zn - Zi Z t - Z 
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Inserting the logarithmic derivative of (|63|l for n ^ m and using (|64(l we find by straightforward compu- 
tations that the condition (|61(l is identically satisfied provided 

Cp n = -M n (66) 

and the functions / and $ satisfy the following functional equations, 

f(x) + f'(y) = 2 f{x) ~J iy \ Z(x) = f(x)f(x) (67) 
x y 

for all x ^ y. The general solution of these functional equations is given by l|51l) . and this concludes our 
proof. Note that at the end we can make the simplifying assignment tp n = M n = 0. 



IV. Alternative formulations 

The strategy to obtain alternative formulations of " goldfish- type" iV-body problems is by now standard 
(and quite old 6 ; for a convenient up-to-date presentation see 2 ). One introduces a monic polynomial 
•4>{z,t) of degree iV in z, the N zeros z n (t) of which evolve according to the equations of motion of 
the iV-body problem under consideration, and then investigates the corresponding evolution of the iV 
coefficients c m (t) of this polynomial. The route we follow to obtain the equations of motions satisfied 
by the coefficients c m (t) - equations that are of course no less solvable than the equations of motion 
satisfied by the zeros z n (t), since the relationship among these quantities, the N zeros z n and the N 
coefficients c m of a polynomial of degree N, is purely algebraic - is via the evolution equation satisfied 
by the polynomial ip(z, t): note that this entails that this evolution equation is itself solvable. Since this 
technique is by now standard, and appropriate formulas to implement it are available (see in particular 
the Appendix in Ref. [2]), we present without further ado the relevant results. 

The evolution equation satisfied by the polynomial ip(z,t) (see Q130 ) the zeros of which evolve according 
to the equations of motion (JJJ reads 



1> tt - 2 (z 2 - a 2 ) 4, tz + 2 [(JV — 2)z — ci] */> t 

+ (z 2 - a 2 ) 2 V zz - 2 [(N - 3) z - ci] (z 2 - a 2 ) ^ 

+ {N (TV -5) z 2 -2 (N-2) Cl z + 2 [2 N a 2 + c 1 - cf + 3 c 2 ] } tp = . (68) 

Notation: here and hereafter subscripted variables denote partial differentiations with respect to them. 
Remark J^.l. This evolution equation, (|68() . contains also certain coefficients c m = c m (t), which are 
obviously (linearly) related to the function ip(z,t), indeed clearly (see l|13|l) 



c m (t) - [(N-m){ 



_! d N - m tp{z,t) 



(69) 



Hence (|68H is in fact a nonlinear functional equation satisfied by the polynomial ip(z,t), and the fact that 
it is indeed satisfied by a polynomial of degree N in z, while not evident, is implied by the way it has 
been obtained. □ 

From this evolution equation one obtains (using if need be the results in [5]) the corresponding system 
of ODEs satisfied by the coefficients c m (t), see (|14f) : and this of course justifies the relevant results about 
the solvability of this system reported in the preceding Section II. 

Exactly the same procedure yields l|21|l from ©, although a more direct route is via the "trick" 
formula (J5J, as indicated in the preceding Section II. Anyway we also display here, for completeness, the 
equation for the polynomial ip(z,t) (see Q22JI1 that provides the bridge connecting these two systems of 
ODEs (deriving this equation is particularly easy using the formulas given in the Appendix in Ref. 
but beware of the slight notational change in the definition of the coefficients c m due to the (i) m factor 
in the right-hand side of 122fl ): 

]>tt-2z (z-i) i> tz + [2 (N-2) z-{2N + l)i-2ici] & 
+z 2 (z — i) 2 ip zz — 2 z (z — i) [N (z — i) — 3 z — i Ci] 4> z 
+ [N (N - 5) z 2 - 2 iV 2 % z - N (N + 1) - 2 (N - 2) i c x z - 2 (N - 1) c x 

+2(ic 1 + c 1 2 -3c 2 ) V> = 0. (70) 
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V. Equilibrium configurations, behavior in their vicinity, Dio- 
phantine relations 

In this section we discuss the equilibrium configurations (namely, the time-independent solutions) of the 
isochronous models (jHJ and l|21[l and the behavior of these models in the vicinity of their equilibria. 
A motivation for focusing on the isochronous models is that they lead to the remarkable Diophantine 
relations reported at the end of Section II, as indicated below. 

Clearly the equilibrium configuration c m (t) — c m , c m (t) = of the system of ODEs (|21J) is character- 
ized by the following system of N algebraic equations: 



- (m + 2) (to - 3) c m +2 + 2 (to - 1) (to + 1 + ci) c m+ i 
+ [-m (m + 1) - 2 (m - 1) ci + 2 c\ - 6 c 2 ] c m = , 
to = 1, N , c = 1 , c_i = c N+ i = c N+2 = . 



(71) 



Likewise the equilibrium configuration z n (t) — z n , z n = 0, of the iV-body problem JBJ) is characterized 
by the following N algebraic equations: 



Z n + i + 



N 

E 



-•n 



= o 



These two configurations are related to each other by the polynomial formula (see J23) 



JV 



$(z) = n (* 



n=l 



n=0 



, AT— m 



CO 



where il>{z) is the "equilibrium" (namely, time-independent) polynomial solution of (|7U|I . 
The general solution of the algebraic problem l|72l) can clearly be broken down as follows: 



E 



= for n = 1, 



m— l,m^n 

z„ = i for n = !/ + 1, .../j, 



(72) 



(73) 



(74a) 
(74b) 



z n =0 tor n = fi + l,...,N , (74c) 

with v and /i nonnegative integers, < i> < fi < N . Of course each of these 3 sectors will be empty if 
the corresponding range of values of n is empty (recall that n — 1, N). 

Remark 5.1. In any equilibrium configuration the labeling of the particles can be freely permuted. To 
write the breakdown (|74|l we identified, without loss of generality, a (somewhat) definite assignment of 
particle labels. □ 

Remark 5.2. Genuine equilibrium configurations of the iV-body problem 10 are characterized by 
the requirement that z n ^ z m if n ^ to: indeed, whenever this condition is violated, the equilibrium 
condition l|72l) becomes ambiguous due to the vanishing of some denominator in the sum, compensated by 
a vanishing of the corresponding numerator or by some other cancellation. Hence a necessary condition 
in order that the configuration associated with the breakdown indicated in (|74fl correspond to a genuine 
equilibrium configuration of the A^-body problem |J?J| is that /i > TV — 1 and v > fi — 1, so that at most one 
of the ZnS vanishes (in which case we assign to it the highest label, zn = 0) and at most one takes the 
value i (in which case we assign to it the highest or next-to-highest label, zn = i, or zn-i = i if zjv = 0). 
But in the following it is convenient to consider all possible equilibrium configurations, including non 
genuine ones, because, as we will see, such configurations, while problematic to deal with in the context of 
the A^-body problem JBJ, correspond to equilibrium configurations c m of the system of ODEs (|21|l which 
are instead perfectly legitimate in the context of this nonlinear harmonic oscillators model. Indeed their 
consideration in such a context yields interesting findings (see below). □ 
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To get more information on the roots z n of (|74a|) we now introduce a monic polynomial (p(z) of degree 
v having the v numbers z n with n — 1, v as its zeros: 

^)=II( z - z "») = E( ! ) m ^ / "'"' ^o = l- (75) 

n— 1 m— 

Note that via this formula we also introduced the v coefficients p m of this polynomial. 

It is now straightforward (and particularly easy using the formulas given in the Appendix in Ref. 
but beware of the slight notational change in the definition of the coefficients c m due to the (i) m factor 
in the right-hand side of J75J) to conclude that this polynomial must then satisfy the following equation, 
implied by (|74a|l : 

z 2 <p" -2 (y - 3) zip' + v {v - 5) p 
= i [zip" -2 (v + ipj) p'] . (76) 

It is now easily seen that, via l|75|l . this ODE 1|76|) yields for the coefficients p m the recurrence relation 

m (m - 5) <p m = (m - v - 1) (m + v + 2 tp x ) p m _ 1 , (77a) 
which must be complemented by the two extremal conditions (see (JE30 

f-i = Vv+x = (77b) 
and by the normalization condition (see l|75|)) 

Vo = 1 • (77c) 

Clearly the two extremal conditions l)77b(l are identically satisfied (for m — respectively m = v + 1), 
while the condition (|77c|) yields (for m = 1) 

entailing the requirement (hereafter assumed to hold) 

v ^ 2 . (79) 
Insertion of l|78|l in (|77a() yields finally the recursion 

m (m - 5) ip m = (m - v - 1) {m + P m -i , (80) 

the solution of which is easily seen to exist only if v < 5. For the remaining cases, 

v = or 1 or 3 or 4 or 5 , (81) 
see (|79|) . we get the following solutions (recall l|77c|l l: 

for v = , p(z) = <Pq = 1 , (82a) 

for = 1 , = Vi = 1; v{ z ) = z + i , (82b) 
for v = 3 , ¥? = 1, p x = -6, p 2 = 14, p 3 = -14, , (82c) 

fori/ = 4, Vm = (-) m Q) 5 m = 0,l,..,4, (82d) 
for z/ = 5 , 99 m = (-) m [ ) , m = 0, 1, 4 , <y9, arbitrary . (82e) 
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Via 174(1 and (|75|l it is clear that the monic polynomial ^p(z) of degree N in z, see (|73|) - which 
identifies as its N zeros z n respectively its N coefficients c m the equilibrium configurations of the models 
© respectively (|2T|) - is given by the formula 

j>(z)=<p(z) {z-iY~ v z N -» . (83) 

ft is thereby seen, via l|82(l . that the coefficients c m are given by the formulas (|28|l (with the arbitrary 
constant c = (p 5 — 1 in l|28c(l . sec (I82e(l 'l. 

Remark 5.3. Clearly the coefficients c m vanish for m > p,, hence they all vanish (except of course 
Co = 1) if p = (this assignment provides indeed a solution of l|71j0 . □ 

Next, let us discuss the behavior of the system (I21|l in the neighborhood of its equilibrium configura- 
tions. To this end we set 

c„(t) = c n + e Pn (t) + 0(e 2 ) , (84) 

with c m the coefficients c m at equilibrium (as determined above) and e a small parameter. We thereby 
obtain in the standard manner the linearized equations of motion 

Pm + 2 ( m - !) « Pm+l - ( 2 m + 1 + 2 5 l) i P m 

- (m + 2) (to - 3) p m+2 + 2 (to - 1) (to + 1 + ci) p m+1 
+ [-m (m + 1) - 2 (m - 1) ex + 2 cf - 6 c 2 ] p m 
+2ic m p 1 + 2 [(to - 1) c m+1 - (to - 1 - 2 ci) c m ] Pi - 6 c m p 2 = , 

m=l,...,N , p = 0, = p N+2 = . (85) 

The general solution of this linear system of ODEs, JHSJ), reads 

N 

Pm® = E K +) r rn Kn) CX P *) + «P ^ eX P foP *)] - ( 86 ) 

n=l 

where the 2 numbers are arbitrary (to be fixed by the initial data) while the 2N numbers pit\ 
respectively the 2N corresponding (^-independent) iV-vectors = ( r^^ n \ r^^ n ^ , are the 

eigenvalues, respectively the eigenvectors, of the (N- vector) generalized eigenvalue equation (|29a|) . This 
implies (|29b|) with the two N x N matrices A and B defined (componentwise) by the formulas 1271) . 
But we know (see Proposition 2.13) that all the nonsingular solutions of the system of nonlinear 
harmonic oscillators (|21|l are completely periodic with period 2 7T, hence the (certainly nonsingular) 
solutions describing the behavior of this system around equilibrium must have the same periodicity 
property, implying that all the eigenvalues yielded by the generalized eigenvalue problem H29a|) must 
be integers. And this entails the validity of Proposition 2.14. 

Remark 5.4- In the special case of the equilibrium configuration c m — (for to — 1, AT, while of 
course Co = 1; see the Remark 5.3), the matrices A and B become triangular and the computation of the 
eigenvalues is then a trivial task, yielding 



P 



W=n + 1, p£->=n.H (87) 



VI. Outlook 

It is remarkable that a research project started with the main purpose to clarify a methodological issue 
- namely, the relationship among two different approaches to the same question: that of identifying 
solvable many-body problems - resulted in the identification of a novel solvable many-body problem. To 
the readers who might imagine - in view of the recent discovery of several such new models, as reviewed in 
Ref. [2] - that this is a relatively trivial task, we suggest to try and find themselves some new model. Our 
educated guess is that such a task is quite challenging. We are nevertheless ourselves hopeful that new 
many-body models exist, and that they might be discovered/manufactured by the techniques described 
in this paper. In any case this possibility remains as a tantalizing prospect, until a way is found to 
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ascertain conclusively that these approaches have exhausted their capability to yield many-body models 
of the kind investigated herein which are both new and interesting (although the second of these two 
qualities involves of course a value judgment). 

Another research direction (perhaps suitable as a PhD project) is towards proving the Diophantine 
conjectures proffered in this paper (see the end of Section II) and in previous ones (see [2] and other 
papers referred to there), as well as obtaining additional findings of this kind (for instance by applying 
techniques analogous to those of Section V to the model Q14JI. taking advantage of the results reported 
in the Appendix). 
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Appendix: Equilibrium configurations of the models ( 1141) and (OQ) 

In Remark 2.9 it was mentioned that there exist additional equilibrium (i. e., time-independent) solutions 
of the model 1|14|) besides l|16(l . In this appendix we firstly list all these equilibrium configurations (as 
obtained via Maple) for N = 2,3,4 and 5. We then outline a technique allowing to obtain all the 
equilibrium configurations for arbitrary N (we of course did check that these findings reproduce, for 
N = 2, 3, 4,5, those obtained via Maple). 

These equilibrium configurations are clearly solutions of the following set of N algebraic equations 
(see ijTi))): 

(to + 2) (to - 3) c m+2 - 2 (to - 1) c\ c m+i 
+2 [to (N + 2 - to) a 2 - c 2 + 3 c 2 ] c m 

-2 (N + l - to) a 2 ci Cm_i + (N + 2 - to) (N+l-m) a 4 c m _ 2 = , 
to = 1, ...,N , Co = 1 , c_i = c N+1 = c N+2 = . (Al) 

The following solutions of this algebraic system have been obtained via Maple. 
For N = 2 

d = , c 2 = -a 2 , (A2a) 

or 

c 2 a 2 

c-2 = — — , ci arbitrary (A2b) 

(the equilibrium configurations (|16|l correspond to the latter one, (|A2b|l , with c\ = ±2 a) . 
For N = 3, 

c? aci 2 ac? 2a 2 ci , . . 

c 2 = — ± — a , c 3 = ±— — , ci arbitrary (A3) 

(the equilibrium configurations (|16|l obtain for ci = ±3 a) . 

For N = 4, 

c i 4 a 2 2 a 2 Ci a 4 , . . . t . 

c 2 = -jj — , c 3 = -a ci , c 4 = Cl arbl trary , (A4a) 



18 



C2 
C4 



2 1/ C] 4> 



2a' 



,2 

3 3 

a 2 c 2 4 a 3 ci 



C3 = ± 



2 a c 2 5 a 2 ci 



a , ci arbitrary 



(the equilibrium configurations (|16fl obtain from (|A4bl) for ci = ±4 a) . 

For TV = 5, 



(-2 
C 5 



-i ± — - - 2a' 



cf , aci 
3 ~ 3 

2 
1 



a 3 c 2 . 2 a 4 ci 



C3 = ± 



ac 2 5a 2 ci 



3 3 
ci arbitrary , 



C4 



a Ci 



(A4b) 



(A5a) 



C2 
C5 



i ±aci 



10 a' 



, a 3 Ci 2 4 8 a 

± 2 a ci ± , ci arbitrary 

3 3 



C3 = ±a — 3 a ci , C4 = a c x =p 5 a ci + 5 a 
5 



(A5b) 



(the equilibrium configurations (|16|l obtain from (|A5bl) for c\ — ±5 a). 

The route we follow to obtain all the solutions of the system (|A1|) is analogous to that followed in 
Section V. The starting point is to introduce a monic polynomial ip(z) of degree N that has the 
numbers c m solutions of (|A1|> as its N coefficients: 



A' 



iV 



, N- 



c = 1 . 



(A6) 



Note the analogy of these formulas with 113(1 . and the fact that we also introduced the N zeros z n of 
this polynomial t/j(z), which clearly provide the equilibrium configuration of the A^-body problem |T]l 
(although not necessarily a genuine equilibrium configuration), hence satisfy the following system of N 
algebraic ODEs: 



(z 



, 2 -a 2 ) 



Zn + 



N 

E 







(A7) 



Our strategy to find all the solutions of the system i|Al|) is to find firstly all the solutions of this 
system, (|A7|) . and then use HA6J1 . 

Clearly the solutions of l|A7|l can be broken down as follows: 



Zn + 



E 



= for n = 1 , . . . , v 



z n = a for n = v + 1, ...,v + /i , 



(A8a) 
(A8b) 



z n — — a for n = v + /! + 1, A^ , 
with the two nonnegative integers v and fi arbitrary except for the constraint 

v + (J, < N 

(implying of course that neither one of these two nonnegative integers can exceed N). 
This assignment, 1A8JI . clearly entails that 



i){z) = (z — a) M [z + a) 



(A8c) 
(A8d) 

(A9a) 
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with <j> v (z) the monic polynomial of degree 

V V 

<!>*(*) =H (*-*»)= f m z v ~ m , f = l, (A9b) 

n— 1 in — 

the zeros of which satisfy the algebraic relations (2BHJ- Hence this polynomial <\> v {z) satisfies the equation 

(z 2 - a 2 ) cf>'l - 2 [{u - 3) z - h) fa + v {v - 5) fa = , (A10) 

as implied by the (by now standard) technique to transform algebraic equations such as IjASaJl into 
differential equations (see for instance the Appendix in Ref. jz]). Here and below primes denote of course 
differentiations with respect to the argument of the function they are appended to. 

Before proceeding to discuss the solution of this equation let us consider the special case with v = 
entailing 4> (z) — 1 (which solves (|A10|) trivially). In this case (|A9a() and l|A8d|) yield 

4>{z) = (z-a) M {z + a) N ~^ , (i = 0,l,...,N . (Alia) 



It is then easily seen from (1A6|) that this entails 

min(/i,m) 



fcmax(0,m+(i-JV) 



This formula provides a set of equilibrium configurations, characterized by the integer /i in the range 
< /i < N; in particular the two solutions corresponding to /i — and to fj, = N are easily seen to yield 
the two solutions 1)16(1 . 

Let us now return to ((AlOjl . assuming hereafter that the integer v is positive, v > (to avoid 
unnecessary notational complications). To solve this equation, 1A10(1 . we set 



<P v (z) = a v X (x) , z = a (x-1) . (A12) 

This formula implies that x( x ) is again a monic polynomial of degree v (although for notational simplicity 
we do not signal this via a subscript v). We also set (in analogy to l|A9b|l 'l 

X (x) = l[(x-x n )=J2x m x v - m , Xo = l, (A13) 

n— 1 m— 

and we then note that this formula, together with (|A9b|l . entails 

m / — f\ 

/■» = ° m £( m _/)*' ( A14a ) 

hence in particular 

h=a(y + xi) ■ (A14b) 
Via this formula and l)A12(l the differential equation l(A10|l now reads 

x (x - 2) x" - 2 [{v - 3) x - 2 v + 3 - X i] x' + v (y - 5) x = , (Al5a) 

entailing, via (|A13|) . the two-term recurrence 

m (m - 5) Xm = 2 { v + 1 - m ) ( 3 - v ~ Xi - m) Xm -i , (A15b) 

implying (for m = and m = v + 1) the extremal conditions X—x = X v +i = (consistently with (|A13|l ). 
For to = 1 the recurrence formula (|A15b|) (together with Xo = 1, sec QA130 ) yields the relation 

(u - 2) xi = -f (i/ - 2) , (A16) 

requiring that the two cases with v = 2 and ^ 2 be treated separately. 
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For v = 2 one easily obtains the solution 

X6 = 1 , xx arbitrary, x 2 = *lM±1) . (A17a) 

Remark A.l. Let us note as a curiosity that the recursion (|A15bJl with v — 2 allows this solution 
(|A17a|) to be extended as follows: 

X 3 = Xi = , X 5 also arbitrary, 

2— 3 15 / m -l + x A 

Xm = — / TT7 , X 5 i m = 5,6, ... . (A17b) 

to [m — 1) (to — 2) \ m — 5 / 

But we are only interested in the solution (|A17a|) with x m — f° r to > 2, entailing that x(x) is a 
polynomial of degree j/ = 2 . □ 
So in this ^ = 2 case we get 

x(o;)^ a ; 2 + x 1 x+ Xl (X 3 1 + 1) (A18) 

hence, via (|A12|I . 

f 2 - a 2 

4> 2 {z) = z 2 + fiz+ , f x arbitrary (A19a) 



where we set, consistently with (|A9bl) . /i = (2 + Xi) 
For v ^ 2 the recursion l|A15b[) with (|A16|) yields 

Xo = 1 . Xi = -v i X 2 = 3 i X 3 = X 4 = . X 5 arbitrary, 

X™ = — ; TT7 o\ [ k Xs ! to = 6, A20 

TO (to — 1) (to — 2) \to — 5/ 

Of course the second line of this equation is only relevant if v > 5, which can only happen if N > 5. 

From these results, via (|A9all with (|A8dl) and l|A12() . we arrive finally, after a bit of trivial algebra, at 
the following two determinations of the polynomial 4>(z) (see l |A61 ): 

i>(z)= (z 2 +cz+ ° 2 ~° 2 ) (z-a) M (z + a) N ^ tl , fi = 0,1, N — 2 , (A21a) 



^(z) = (z-a) M (js + a)* - " 
(-2)' 



v a v (y — 1) a' 



+c 



v-5 
£=0 



z + a 3 (z + a) 



i/-5 



a 



+ 5) (* + 4) (£ + 3) V * / V^ + « 



= 5, 6, N , ii = 0, 1, N - v 



(A21b) 



The first, l|A21ali . of these two formulas is applicable for N > 2, having been obtained from the previous 
results corresponding to v — 2, with c = fi an arbitrary number. It includes the results obtainable from 
the cases with v = 1 and v = 3; likewise, the result corresponding to v = 4 has not be reported, as it is 
encompassed by the result (|Alla|) . The second, <|A21b|) . of these two formulas is of course only applicable 
provided N > 5, having being obtained from l|A20fl for v > 5, with c = X5 / 60 an arbitrary number. 
Together with IjAllajl these two formulas determine all the equilibrium configurations c m of the system 
(114(1 (and as well all the equilibrium configurations - not necessarily genuine - of the A-body problem 
0) - up to a final (trivial but tedious) step, to be performed using l|A6(l (as done above to obtain (|A1 lbfl 
from HAllaj l). which we leave as a task for the diligent reader. 
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